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ABSTRACT 



We present the first results of a survey of random fields with the slitless G141 
(Ac = 1.5yU, AA = 0.8/i) grism on NICMOS. Approximately 64 square arcminutes have been 
observed at intermediate and high galactic latitudes. The 3o" limiting line and continuum 
fluxes in each field vary from 7.5 x 10~^^ to 1 x 10^^^ erg cm~^ sec~^ and from H = 20 
to 22, respectively. Our median and area weighted 3cr limiting line fluxes within a 4 pixel 
aperture are nearly identical at 4.1 x 10~^^ erg cm~^ sec~^ and are 60% deeper than the 
deepest narrow-band imaging surveys from the ground. We have identified 33 emission-line 
objects and derive their observed wavelengths, fluxes and equivalent widths. We argue that 
the most likely line identification is Ha and that the redshift range probed is from 0.75 to 
1.9. The 2cr rest-frame equivalent width limits range from 9A to 130A with an average of 
40A. The survey probes an effective co-moving volume of 10^ h^^ Mpc^ for go = 0.5. Our 
derived co-moving number density of emission line galaxies in the range 0.7 < z < 1.9 is 
3.3 X 10~^ /150 Mpc~^, very similar to that of the bright Lyman break objects at 2; ~ 3. The 
objects with detected emission-lines have a median F160W magnitude of 20.4 (Vega scale) 
and a median Ha luminosity of 2.7 x 10*^^ erg sec~^. The implied star formation rates range 
from 1 to 324 Mq yr~^, with an average [NII]6583,6548 corrected rate of 21 Mq yr"-*^ for 
Ho = 50km/s/Mpc and go = 0.5 (34 Mq yr~^ for go = 0.1). 



-3- 
1. Introduction 

Spectroscopic surveys of galaxies at 2; < 1 have reached a reasonable level of maturity 
in recent years (e.g. CFRS, Lilly et al. 1996; Le Fevre et al. 1996; CNOC, Yee et al. 1996). 
The use of specialized selection techniques has led to a similar or greater gains at redshifts 
of approximately 3 and larger (Steidel et al. 1996; Dickinson 1997; Hu, Cowie, & McMahon 
1998), with a few examples being identified to 2; ~ 5 (Dey et al. 1998; Weymann et al. 
1998; Spinrad et al. 1999). However, to build a coherent picture of galaxy formation and 
evolution we are hindered by two problems. The first is that we know little of properties 
of normal galaxies in the region between 1 < z < 2. Secondly, most high redshift galaxies 
now known are selected at rest-frame UV wavelengths, where the degree of extinction is 
uncertain (e.g. Heckman et al. 1998). The detection of heavily dust-reddened galaxies at 
2; ~ 2 — 3 in the sub-millimeter wavelength (Small, Ivison & Blain 1997) suggests that high 
redshift galaxies selected by UV colors represent only part of the population at these early 
epochs. At present, the estimates of the global star formation rates are uncertain as they 
are based primarily on UV selected samples (Pettini et al. 1997; Meurer et al. 1997). 

One approach to dealing with both the scarcity of spectral signatures in the visible 
bandpass and the uncertainties arising from extinction is to observe in the near-IR. Until 
recently the available detectors and spectrometers have made spectroscopic surveys in 
the near-IR impractical and so most programs have been based on narrow band imaging 
(AA ~ 1%) to search for strong emission-line objects, primarily in the K window where Ha 
can be detected for 1.9 < z < 2.6. Two narrow band near infrared imaging surveys, covering 
several hundred square arcminutes, were carried out with the Calar Alto 3.5m telescope 
(Thompson, Mannucci & Beckwith 1996; Beckwith et al. 1998; Mannucci et al. 1998). The 
first Calar Alto survey targeted fields containing known quasars or radio galaxies and the 
second the redshifts of damped Lya or metal line absorption systems. The first survey 
discovered one emission line galaxy aX z = 2.43 over an area of 276 square arcminutes to a 
3a fiux limit of 3.4 x 10~^^ erg cm~^ s~^, while the second survey found 18 emission line 
galaxies (16 objects with detected Ha and 2 with [OH] 3727) over an area of 163 square 
arcminutes to a 3a fiux limit of 2 x 10^^^ erg cm^^s^^. To date the deepest narrow band 
near-IR imaging survey is that carried out by Teplitz, Malkan & McLean (1998) with NIRC 
on the Keck 10m telescope. They found 13 emission line galaxies in an 11 square arcminute 
area to a 3a fiux limit of 7 x 10^^^ erg cm^^s^^. The target fields of this survey were also 



-4- 



selected by matching the redshifts of quasar emission and metal absorption hnes to their 
narrow band filters. The space densities of strong Ha emitting galaxies inferred from these 
surveys span more than 2 orders of magnitude. The large disparities in depths between 
the Calar Alto and Keck surveys and the strong biases introduced by selecting redshifts 
matched to known absorption-line clouds make comparisons between the surveys and with 
other measures of the average star formation rate non-trivial. 

NICMOS offered a unique opportunity to address this problem with multi-object 
spectroscopy. The wide field camera on board NICMOS (Thompson et al. 1998) is 
equipped with three grisms that allow slitless spectroscopy over 0.75 square arcminutes per 
exposure. The extremely low background at wavelengths below l.dfi allows for a sensitive 
low resolution spectroscopic survey for Ha and other emission-lines through volumes that 
are an order of magnitude larger than those available to narrow-band imaging surveys due 
to the large redshift range sampled by the grism. 

Throughout this paper, Hq = 50km/s/Mpc and go = 0.5 are adopted unless explicitly 
indicated otherwise. 



2. Observations 

All of the data presented here were obtained with Camera 3 on NICMOS using the 
G141 grism and the F160W and FllOW broad-band filters. The G141 spectral element 
covers the wavelength range from l.l/x to 1.9yU with a mean dispersion of 8.1 x 10~^/i 
per pixel. The nominal resolving power of the grism is 200. The actual resolution that 
is achieved is a function of the apparent image size and the camera 3 PSF. With the 
PSF present in the parallel observations, the realized resolution is rarely better than R 
~ 150. Small variations in the PSF due to changes in the optical telescope assembly (e.g. 
breathing) and longer term changes in the internal structure of NICMOS introduce small 
variations in the maximum achievable resolution. As discussed in detail below many of 
the objects detected in our survey are spatially resolved and the resolution of their G141 
spectra is determined by the image size. 
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2.1. Parallel observing mode and scheduling 

All of the observations for this investigation were obtained with HST operating in the 
parallel mode. While this allowed us to collect far more data than would be possible in 
a single primary program, it limited our ability to plan and execute the observations in 
a manner that optimized the science return. Useful data for this program first became 
available in October 1997 when the default position of the pupil alignment mirror was 
moved from the camera 2 optimized position to the best reachable position for camera 
3. After that time camera 3 parallels were scheduled during much of the time for which 
WFPC2 or STIS were the primary instruments. Most of the data in this investigation was 
obtained as part of a public service parallel program that was implemented and managed by 
the Space Telescope Science Institute. The observing algorithm was quite simple. In each 
full orbit one of the following exposure sequences was selected: F160W and FllOW imaging 
with exposure times of either 512 or 382 seconds, F160W (2 x 256 sec) and G141 (2 x 896 
sec), or FllOW (2 x 191 seconds) and G096 (2 x 896 sec). The selection of the exposure 
sequence for any given orbit was nearly random, but was weighted in favor of the purely 
imaging sequence for much of the October 1997 to March 1998 period. After that time the 
balance was shifted towards spectroscopy with the grisms. All of these observations were 
planned and implemented by the STScI staff, with the guidance of Dr. John MacKenty. 

The array was read using the STEP64 read pattern. A 2 dither step was executed 
between each pair of direct or grism exposures. This dither step was made in the detector 
reference frame and was along array rows. The dither step and orientation were fixed, 
resulting in only two independent positions regardless of the number of exposures for any 
given target. Small moves executed as part of the primary observing plan introduced 
additional offsets for some of the longer pointings. 

Approximately 10% of our data were obtained as part of a crafted guest observer 
(GO) parallel program that was implemented for a portion of the February 1998 to May 
1998 period. The basic approach was similar to the public parallel program, except that 
the dithering strategy was more complex. For each pointing we began with a number of 
F160W images with exposure times selected to equal 15% of the expected total duration of 
that particular pointing. Dither moves were executed between each exposure, but with step 
sizes and orientations that were non-redundant. The remainder of the orbit and subsequent 
orbits were devoted to G141 exposures with a dither pattern matched to that of the direct 
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images. For a given exposure time the noise levels in these crafted observations were ~ 20 
to 30% lower than those in the public parallels. 



2.2. Data processing and analysis 



2.2.1. Two dimensional spectra 

We began our treatment of the data using the calibrated output from the STScI 
data processing pipeline. The pipeline reduction consists of dark current removal on a 
read-by-read basis, linearization, bad pixel masking, and cosmic ray rejection. At this stage 
the data remain in the format of two dimensional, single frame images for each exposure. 
The principal elements of the reduction process that are accomplished by the pipeline to 
this point are the first order bias and dark subtraction, linearization and CR rejection. Our 
subsequent processing improves on the bias and dark subtraction, bad pixel masking and 
more complete rejection of particle events. The calibrated grism exposures have not been 
corrected for the flat field response and have not had the sky background removed. The 
flat-field response is a function of wavelength for each pixel and so depends on the location 
of each object, and will not be the same for the object and the sky. The flat-fleld correction 
is applied during the 1-D spectral extraction. 

We constructed a median sky frame for the G141 grism by combining the calibrated 
exposures for approximately 200 frames, covering 80 independent positions on the sky. We 
rejected observations made at low galactic latitudes as well as those centered on bright 
objects or crowded flelds. This produced a high signal-to-noise uncorrected sky frame 
for use in the sky subtraction for each exposure. The observed background rates vary 
signiflcantly from frame-to-frame depending on the ecliptic latitude and longitude. We 
subtracted a scaled version of the median sky from each calibrated 2-D spectral frame. 
The optimal scaling for the sky was found by minimizing the RMS noise in two object-free 
regions of each frame. Very low-spatial frequency structure remained in the background of 
many frames, primarily due to the shading effect in the NICMOS arrays. We removed a low 
order surface from each frame. This surface was derived from a second order polynomial flt 
to an intermediate image constructed from the median of each column. This surface fltting 
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was carried out separately in each quadrant of the detector, allowing us to remove small 
quadrant-to-quadrant bias variations. The individual sky subtracted and shading corrected 
frames were then shifted into registration, using offsets derived from the accompanying 
direct images. The registered 2D spectra were then averaged with a bad pixel mask and 3a 
clipping applied to remove hot pixels and residual cosmic rays. 

The direct images that are used to calibrated the extracted one dimensional spectra 
were processed and combined using the same techniques as above, but without the sky 
subtraction and low order surface fitting correction steps. 



2.2.2. 1-D spectral extraction and calibration 

The extraction and calibration of the one-dimensional spectra was carried out with the 
NICMOSlook software developed at ST-ECF (Freudling et al. 1998). The extraction was 
done in the following war. First, the position and size of the objects previously identified 
as spectral line objects were interactively marked on the image display. The extraction 
region of the spectrum was then computed using transformations derived as part of the 
calibration program. The extracted spectra were corrected for the wavelength dependent 
flat-field response on a pixel-by-pixel basis. For each pixel, the quantum efficiency for the 
extracted wavelength was estimated by interpolation from a series of narrow-band flat field 
exposures. Fluxes were calibrated by applying throughput curves previously derived from 
the observation of fiux calibrators. The final one dimensional spectra were obtained by 
adding the fiux of all pixel which correspond to each wavelength bin. For some objects 
close to one of the edges of the grism images, only partial spectra could be extracted. Since 
there is a significant displacement between the location of objects on the direct images and 
their dispersed spectrum on the grism image, spectra for a few objects that fall outside the 
area covered by the direct images were extracted. In that case, the wavelength scale was 
estimated from the cut-off of the throughput curves at high and low wavelengths which can 
easily be identified. For such objects, the wavelength scale is uncertain by about O.Ol/x. 

The rms accuracy of the wavelength calibration for G141 is about 2 x 10"'^ /i. Repeated 
calibration exposures have shown no changes larger than 3% in the wavelength zero point 
or 1.5% in the dispersion. The fiux calibration is based on observations of both white 
dwarfs and solar analogs during the NICMOS camera 3 campaigns. The accuracy of the 



spectrophotometry is limited by variations in the quantum efficiency within individual 
pixels. We estimate that the flux scale for an individual spectrum is accurate to about 10%. 

Spectra were also extracted in a parallel effort by members of the team at the NASA 
Goddard Space Flight Center (GSFC) using software independent of the NICMOS-Look 
routines. The GSFC software followed a similar strategy to NICMOS-Look. Spectral 
positions on the chip were determined by the known offsets from the object location in the 
direct image. Extractions were performed on individual (unregistered) frames ffist, with 
registration left to the end. Each spectrum was extracted in a weighted box, similar to the 
procedure for the NICMOS-Look point source extraction. The wavelength solution was 
confirmed by independent analysis of the same calibration data. Flat-fielding was performed 
using the same data-cube of narrow-band observations, interpolated to the wavelengths of 
the extracted pixels, and the flux calibration was performed based on independent analysis 
of the white dwarf observations. Finally, spectra were combined in wavelength space. 
Emission-line wavelengths derived from spectra extracted with the GSFC software were 
found to agree with the NICMOSLook results to within one pixel (O.OOS^m). Line fluxes 
and equivalent widths were measured using Gaussian fits to the line and polynomial fitting 
of the underlying continuum. The overall results of the extraction and analysis using the 
GSFC and ECF/NICMOSLook software were similar. 



3. Results 



3.1. The properties of the data and the identification of emission-lines 

There are several aspects of the data which negatively impact our ability to identify 
genuine spectral features. Since we are observing in the slitless mode, each object produces 
zero, first and second order spectra. The first order spectra contain the useful data, the 
zero and second orders are sources of confusion, the zero order images as they can can 
be mistaken for emission features, particularly when they fall on the first order continua 
of other objects. The zero order images are slightly dispersed and for point sources often 
appear bimodal. The displacement between the zero order images and the start of the 
first order spectrum is 27 , or 135 pixels, a substantial fraction of the array size. For this 
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reason the zero order images are confined to roughly half of the detector area and in most 
cases they can be identified by matching them with either the first order spectra or with 
an object in the direct image. There is a small portion of the detector in which zero order 
images can appear without either first order spectra or images in the direct frame. In 
addition to the confusion caused by the zero order images and second order spectra, there 
are artifacts associated with hysteresis in the HgCdTe detectors. Bright stars imaged with 
the broad-band filters just before a grism observation will produce a low level after-image 
that is displaced from the first order spectrum. Since the data are nearly always taken in 
pairs of direct and grism exposures, these persistent images appear in multiple exposures 
and can take on the appearance of isolated emission-lines. Strong cosmic rays events have 
similar properties, but the two position dithering sequence aids in their identification and 
removal. Particularly strong particle events effect more than one pixel and can produce 
amorphous clouds of electrons that decay over time scales of an orbit. During the South 
Atlantic Anomaly (SAA) passage the particle event rate is quite large and most NICMOS 
images or spectra obtained at this time are not usable. The hysteresis in the detectors often 
renders subsequent exposures useless as well. For a significant number of our short pointings 
one of the two 896 second grism exposures has greatly enhanced noise from residual CR 
events associated with the SAA passage. 

Given the wide variety of point-like features that appear in the two dimensional grism 
images we have made our emission-line identifications on the basis of visual inspection. The 
eye can detect compact features to quite low signal to noise levels. To convince ourselves 
that any possible emission line was real we proceed to eliminate all possible false signals, 
zero order images, persistent images, CRs etc. To facilitate the detection of incompletely 
rejected cosmic rays and other noise features we split each observation into two frames 
of equal length and compared the position of possible emission features. Those features 
that appeared in only one of the data sets, or did not move with the dither pattern, were 
rejected. In Figure 1 we show two pairs of F160W/G141 frames and mark examples of each 
of the sources of confusion as well as real emission features. 



-10- 



3.2. Areal coverage and depth 

During the period from October 1997 to October 1998, 3500 camera 3 exposures were 
taken as part of the pubhc and our crafted GO paraUel programs. These contain a mix of 
direct exposures with the F160W and FllOW fihers and spectroscopic exposures with the 
G141 and G096 grisms. A total of 456 pairs of F160W/G141 frames were taken during this 
time. These were distributed in 107 independent pointings with exposure times ranging 
from 1790 seconds to 21,480 seconds. The G141 grism spectra can provide useful data for 
our purposes to quite low galactic latitudes. Rather than adopt a strict latitude cutoff we 
chose to reduce all of the data and reject those with high stellar densities. After some 
experience we rejected all fields with stellar densities greater than 50 per square arcminute 
above a lOo" threshold in the F160W direct image. We are left with 85 pointings, covering 
64 square arcminutes. 

The final depth achieved varied considerably from field to field. For the 1790 second 
exposures the noise level varied by a factor 3. Some of this variation was due to genuine 
background variations, but much of it was due to various non-Poisson sources of noise 
(e.g. SAA passage). For longer pointings the depth increased roughly as the square root 
of the exposure time. We define our limiting depths as 3cr within a 4 pixel (0.16 square 
arcseconds) aperture. This area is a reasonable representation of the area of the camera 3 
PSF at the focus achieved outside of the campaign mode. In Figure 2 we show a histogram 
of the 3o" 4-pixel limits for the G141 data. 



3.3. Detected Emission Features 

In Table 1 we list all of the fields for which we have reduced and examined the G141 
grism data. For each pointing we give the J2000 lAU designation, the galactic latitude and 
the exposure times for the F160W and G141 filters. 

In Figures la and lb we show two typical pairs of direct and 2-D grism images. These 
two grism images show many of the features discussed in the previous section (zero order 
images, first and second order spectra, after images) as well as one emission-line object each. 
After careful examination of all of the spectra we assembled a catalog of 33 emission-line 
detections. Table 2 gives the field names, object designation, J2000 coordinates, observed 



- 11 - 



L 0055+8518 



Fie0W 






33t3+8535 



F160W 











Gltl 






Z 


- 


-- 


^^"^ 


^--- 


z 


CR 


i 




1 EMISSION LINE 






T=- 


.—.. '-^— 


1 




Z 




- 


fli- -'^-^^'^-^■■: 


z 




• 


^^H 






* 






z 


CONFUSION 



£ "» 1=' ORDER 


GHl 




- - "• 








1 EMISSION LINE 

1 


2 " ORDER 


1=T ORDER 



Fig. 1. — Two typical pairs of direct and 2-D grism images. Each field is 51 x 51 and 
oriented in arbitrary directions. In the G141 frames we have marked zero order images (Z), 
residual after images (AI), persistent cosmic rays (CR), first and second order spectra and 
our identified emission- lines. The region marked as "confusion" in the upper panel identifies 
the portion of the detector within which zero-order images neither have corresponding first- 
order spectra within the G141 field of view or direct images in the F160W field. We have 
marked the objects producing the emission lines in the F160W images. 
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wavelengths, redshifts, line fluxes, significance levels, observed equivalent widths and 
equivalent width limits {2a) for these 33 emission line galaxies. We set a rejection threshold 
at 3cr, the derived significance levels of our cataloged emission-lines are given in Table 
2. These 33 emission-lines are distributed in 22 fields, 6 of which contain more than one 
emission-line object. We do not convincingly detect more than one line in any object. 
All of lines are quite strong, in terms of equivalent width. This is expected given the low 
spectral resolution of the grisms; weak lines against bright continua are not easily detected. 
The observed equivalent width limits reflect the signal-to-noise ratios in the continua as 
well as the galaxy sizes. The 2 a observed equivalent width limits are calculated as follows: 
W/jm = 2 X (^-) X vN, with s the scale in A per pixel, snr the signal-to-noise ratio per 
pixel of the continuum at the emission line and N the number of pixels under each line. 
For lines that are not velocity resolved, the apparent line width is the convolution of the 
instrumental spatial resolution and the intrinsic size of the emission line region. For our 
emission line galaxies, we take N to be 3 pixels. The 2a rest-frame equivalent width limits 
range from 9 A to 130 A with an average of 40 A. In Figure 3, we show direct images in the 
F160W band (~ H band) and 2D grism spectra for all 33 emission line galaxies. Figure 4 
presents the extracted one dimensional spectra of all 33 emission line galaxies. 
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Fig. 3. — A mosaic of F160W band direct images and 2D grism spectra of 33 emission line 
galaxies found in our survey. The emission-lines in each object are marked with vertical 
ticks. For a few of the objects the continuum images fall near the edge or just off the edge 
of the F160W field of view. For these objects (0917+8142C, 0931-0449A) we do not have 
F160W images. 
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Fig. 4. — The extracted 1-D spectra of the emission line galaxies. For each object we 
have marked our candidate emission-lines with arrows below the line. We plot the entire 
range of the G141 grism for each spectrum even though parts of the spectrum have fallen 
beyond the field of view of the detector for several objects. The redshifts, assuming an Ha 
identification for the line are given for each object. There are deficiencies in the background 
removal for several objects that lead to significant uncertainties in the continuum level (e.g. 
1237+6215A) 
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We believe that most, if not all, of the emission lines that we have detected are either 
Ha or an unresolved blend of Ha+[NII]6583/6548. In Table 3 we list the line luminosities, 
star formation rates derived from the Kennicutt (1984; 1992) relation, the apparent 
H-band magnitudes (on the Vega scale) and J— H colors for the fields where the FllOW 
observations were available. The identification of the emission-lines as Ha is based on a 
number of arguments. First, the most plausible identifications are Ha+[NH]6583/6548, 
[OHI]5007,4959, H/34861 and [OH]3727. At the peak of the G141 transmission these 
correspond to redshifts of 1.28, 2.0, 2.1 and 3.0, respectively. Studies of nearby galaxies show 
that Ha+[NH] is usually the strongest fine, the second being [OII]A3727 (Kennicutt 1992). 
The equivalent width of H/3 scales with the equivalent width of Ha as W(H/3)= 0.15W(Ha) 
- 4. Even in the bluest Magellanic irregulars and starburst galaxies, W(H/9) is only ~ SA. 
With a resolution of 5000 km sec~^, it would be quite difficult to detect such lines in our 
data. Redshifts of ~ 3 are highly unlikely given the median F160W magnitude of 20.5, 
as shown in Figure 5. The z ~ 3 Lyman break galaxies have typical F160W magnitudes 
that are 2 magnitudes fainter (H ~ 22.5) (Giavalisco et al. 1998). The space density of our 
emission-line objects is sufficiently large that they are unlikely to represent the extreme tail 
of the luminosity distribution of the z ~ 3 Lyman break objects. If not Ha, it is more likely 
that the fines that we see are [OIII]5007,4959 rather than [OII]3727. The [OIII]5007,4959 
lines could be marginally resolved in our spectra, but only if the objects are spatially 
compact. At 2; ~ 1.5, if we take an R — H color of 3, the median H magnitude of our objects 
corresponds to M/j = —21.7 for Hq = 50 km/s/Mpc, qo = 0.5, quite close to MJj determined 
from the Las Campanas Redshift Survey (Lin et al. 1996). In Figure 6, the left panel 
shows the distribution of Ha line-luminosity; the median line-luminosity is 2.5 x 10"^^ erg 
s~\ close to the value of L^^, derived in local surveys (e.g. Gallego et al. 1995, Tresse & 
Maddox 1998). If the emission lines are actually [OIII]5007,4959, the implied rest-frame 
H-band luminosities and emission-line luminosities would be a factor of 2.5—3 larger and the 
likelihood of such objects appearing in fiux limited samples is suppressed by the exponential 
cutoff in the luminosity function. Strong [OIII]5007,4959 emission is more often identified 
with AGN than star-forming galaxies, making it even less likely that [OIII] emitters at 2; ~ 2 
dominate our detections (see section 4). Furthermore, in the redshift range between 1.6 and 
1.8, both [OIII]5007,4959 and [OII]3727 are visible in the G141 passband. Lines with longer 
rest-wavelengths are possible, but unlikely identifications. HeII10830 implies redshifts that 
are implausibly low (~ 0.38) and for most redshifts both of the [SIII]9069,9545 lines would 
be within the passband. Teplitz et al. (1999) have confirmed the identification of Ha in 
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two of our objects (J0055+8518a, z = 0.76, J0622-0018a, z = 1.12) by detecting [OII]3727 
at the appropriate redshifts in spectra taken with LRIS on the Keck lOm telescope. 

Ha is blended with [Nil] 6583,6548 in our spectra. To correct for the contribution from 
the [Nil] doublet, we adopted [NII]6583/Ha = 0.3, the average ratio for all types of galaxies 
in the Gallego et al. (1997) objective-prism survey. Taking [NII]6583/[NII]6548= 3, we 
have /(Ha) = 0.71 x /((Ha + N[II])). Thus the correction to the estimated star formation 
rates in Table 3 is 29%. We have not applied any extinction corrections to either the Ha 
luminosities or the star formation rates. 



3.4. Notes on Individual Fields 

In Figure 3 we present cut-out images and 2-D spectra of all of our emission-line 
candidates. Flux and wavelength calibrated 1-D spectra are shown for the same objects in 
Figure 4. We designate each object by its lAU J2000 field designation and a letter. There 
are a few fields that have noteworthy features and we detail these below. 

0040+8505 A: This is one of the weaker lines, in terms of Wa, in our sample. There are 
only two exposures of this field, but the emission feature is seen in both frames. 

0613+4752 A: The emission-line in this object is more compact than the continuum. As in 
the case of 0040+8505 there are only two exposures for this field, but the narrow emission 
feature is seen in both frames and moves with the dither pattern. 

0622-0018 A: This low redshift object is detected in a very low latitude field (b= —6°). 
The emission-line is spatially extended and appears to follow the continuum distribution as 
revealed in the F160W image. 

0738+0507 A: This is also a low latitude field and the object in question is the brightest 
in our sample. The large equivalent width and compact morphology of the emission- line 
suggests that this object might be an AGN. 

0738+0507 B: The putative emission-line is quite faint and diffuse. There are 6 exposures 
of this field and the emission feature appears with equal strength in the two independent 
subsets of the data. 

0741+6515 A: The emission-line in this object is spatially resolved and its morphology is 
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Fig. 5. — The distribution of apparent H magnitude of emission line galaxies. The 
magnitudes are on the Johnson system. 
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Fig. 6. — The distributions of redshifts and Ha line luminosities of the emission line galaxies. 
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similar to that of the continuum shown in the F160W image. 

0741+6515 B: The emission-hne in this object is at the extreme red end of the continuum 
spectrum and is spatially offset from the continuum. The emission-feature is located in a 
portion of the field such that it cannot be a zero-order artifact. 

0917+8142 A: This emission feature is also at the extreme red end of the continuum. At 
the location of this feature there is the possibility of zero order images. If this feature were 
a zero order image, the object would be just off the edge of the F160W direct image. The 
first order continuum of a zero order image at this location would be partially on the G141 
frame, but there is a confusing signal from a bright galaxy at this location. The spatial 
coincidence between the continuum and emission features led us to classify this as a real 
emission-line, although this object should be considered less secure than the others. 

0917+8142 C: The putative emission-line in this object is at the edge of the G141 field of 
view. The feature is only present in one of the two dithered frames, it was moved off the 
detector by the 2 dither motion. The object is also located just outside the F160W field in 
both dither positions. The great strength of this line makes it a candidate AGN. 

1120+2323 A,B: The two emission-lines in this field are extremely compact. Both emission 
features are seen in the two dither positions. For reasons that are not entirely clear, the 
focus in this particular set of exposures was unusually good. The difference may have arisen 
from breathing of the telescope optical assembly, or from motions within the NICMOS 
dewar. 

1143-8036 A,B: This pair of emission-line galaxies appear to have nearly identical redshifts 
in the 2-D G141 frame. The slightly different locations in the F160W image result in 
different wavelength scales, and hence different redshifts. The close proximity of the 
two images make accurate centroiding of the F160W images difficult and this introduces 
uncertainties into the wavelength calibrations. 



4. Discussion 

This survey provides a unique window on galaxies at intermediate redshifts. In so 
doing it provides us with a view of the bright end of the Ha luminosity function averaged 



-22- 



over a cosmologically significant volume. From these data we will derive an Ha luminosity 
function and a measure of the volume averaged star formation rate that is fairly insensitive 
to the effects of extinction. These analyses require careful treatment of the selection 
efficiencies and biases inherent in a slitless low-resolution survey and are presented in Yan 
et al. (1999). In the present work we restrict ourselves to a first-order census of bright 
emission-line objects and a brief comparison with other surveys in this redshift range. 

There are a number of near-IR emission-line imaging surveys that are either in progress 
or have been recently completed. In Table 4 we compare the characteristics of our survey 
with several ground-based programs. Clean comparisons are difficult since our survey is of a 
different character and spans a very wide range of exposure times from field to field. In our 
deepest pointings we achieve fiux limits that are substantially deeper than the narrow-band 
imaging survey being carried out by Teplitz et al. (1998) with the Keck Telescope. Our 
median and area weighted sensitivities are both 4.1 x 10~^^ erg cm~^ sec^^ {3a, 4 pixels), 
~ 60% deeper than the area- weighted detection limit reported by Teplitz et al. (1998). The 
two principal differences between our survey and those listed in Table 4 are the volumes 
probed and the choice of sight-lines. Because our survey is spectroscopic we are able to 
probe large columns in the z direction. The comparable areal coverage of our survey then 
leads to total co-moving volumes that are one to two orders of magnitude larger than those 
surveyed by the narrow-band imaging programs. 

Inspection of Table 4 reveals that most of the narrow-band imaging surveys are 
targeted at lines of sight and redshifts that are a priori known to contain either luminous 
galaxies or absorption line systems. The motivation for this can be clearly seen in the 
results of the first Calar Alto survey (Thompson et al. 1996). Surveys along sight-lines 
without absorbers, or in completely random areas, yield few strong line emitters in the 
volumes that can be probed by imaging with 1% filters. Thompson et al. detected a single 
object in 30,000 h^^Mpc^, implying a space density of 3 x 10^^ /iggMpc"^. The choice of 
AGN and/or absorber redshifts is effective in producing detections, but it clearly impacts 
the inferred space density of line emitters. Teplitz et al. (1998), targeting a mix of QSO 
and absorber redshifts, report the highest space density, at approximately 10~^ /ifgMpc"^. 
This is quite different from the results of the Thompson et al. (1996) survey, even when 
allowing for the very different depths of the two surveys. The second Calar Alto survey 
(Mannucci et al. 1998) demonstrates the effect of targeting metal line absorption systems, 
while Bechtold et al. (1997) have achieved similar results using damped Lya redshifts. 
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Our Ha survey probes an effective co-moving volume of 1.0 x 10^ K^q Mpc^ 
(2.1 X 10^ h^Q Mpc^ for qo = 0.1). Tlius tlie co-moving number density of emission-line 
galaxies averaged over our volume is 3.3 x 10~^ h^^ Mpc^. In comparison, the present-day 
co-moving number density of galaxies brighter than L* in the B band is 7 x 10~^ /i|q Mpc^^ 
(Ellis et al. 1996) and the bright Lyman break galaxies at 2; ~ 3 reported by Steidel et al. 
(1996) have a space density of 3.6 x 10^^ /ijq Mpc^^. The duty cycle of strong Ha emission 
in galaxies with bursting or variable star formation histories can be quite low, given the 
sensitivity of strong line emission to only the most recent star formation. Thus it is not 
unreasonable that our space density is below that of present day L* galaxies and comparable 
to that of the UV continuum selected galaxies at high redshift. The low resolution of the 
NICMOS grisms imprints a strong selection effect against emission-lines with low equivalent 
widths. Given the incompleteness at low equivalent widths and the rather bright H-band 
magnitudes of the detected emission-line objects it seems likely that we are sampling star 
formation in a population that is closely related to the majority population of massive 
galaxies. 

The average star formation rate in our sample is 30 M© yr~^ (40 Mq yr~^ for (go = 0.1). 
After correcting for [NH] contamination, this rate is reduced to 21 Mq yr^^. We have not 
applied any extinction correction to the calculations of Ha luminosities and star formation 
rates. While our sample is not magnitude limited, the distribution of apparent magnitudes 
is sufficiently narrow to allow a reasonable comparison with magnitude limited samples. As 
noted above, our median F160W magnitude of 20.5 is quite close to M|j as derived locally. 
It is interesting to note that the average star formation rate in the Lyman break objects, 
before corrections for extinction are ~ 8 M© yr~^, and that the extinction corrected rates 
are a factor of 3 — 7 higher. Recent near-infrared spectroscopy of 8 CFRS galaxies at z ~ 1 
by Glazebrook et al. (1998) reached an Ha luminosity of 2 x lO^^erg/s/cm^ with a star 
formation rate of 22Mq yr^^, comparable to what we find, although our survey covers a 
wider range of Ha luminosities and star formation rates. 

One source of uncertainty in using our data to derive global star formations rates is the 
unknown degree of contamination from AGNs. The fraction of AGNs in our sample could, 
in principle, be calculated using the AGN luminosity function in the redshift range from 0.7 
to 1.9. However, at z= 1.5 an H magnitude of 20.5 corresponds roughly to an absolute B 
magnitude of -22.5 adopting a B— H color similar that for QSOs observed by Neugebauer et 
al. (1986). At this faint magnitude, the quasar luminosity function is completely unknown. 
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The most recent quasar luminosity function at z ~ 0.5 — 2.0 reaches only to M^ ~ —25 
(Hewett, Foltz & Chaffee 1993). Furthermore, our knowledge of the Seyfert luminosity 
function is entirely based on local samples and the degree of luminosity evolution in the 
Seyfert galaxy populations as a function of redshift is unknown. A crude estimate of the 
fraction of AGNs in faint redshift surveys can be derived from spectrophotometry of the 
CFRS sample. Using various emission line ratios, Hammer et al. (1997) and Tresse et al. 
(1996) estimated that about 8—17% of their galaxies at ^ < 0.7 are active galaxies, either 
Seyfert 2 or LINERs. A contamination rate at this level is consistent with our classification 
of two of our objects as candidate AGNs on the basis of their large equivalent widths 
and compact morphologies. If we exclude the two AGN candidates, the average SFR is 
28.5 M0/yr and 20 Mq/jt if we correct the [Nil] 6583,6548A contamination. The impact of 
line misidentifications is comparable to that of the AGNs, as the most likely contaminant 
are [OIII] 5007,4959 emission from compact non-stellar sources. 



5. Summary 

We have presented the basic data derived from a spectroscopic survey of field galaxies 
in the near infrared. The G141 grism on NICMOS provides a unique opportunity to survey 
large volumes to faint flux levels in Ha and other emission-lines. Our survey of 64 square 
arcminutes at high and intermediate galactic latitudes resulted in the detection of 33 
emission-line objects and we believe that most or all of these are Ha at redshifts between 
0.75 and 1.9. The line luminosities, equivalent widths and continuum magnitudes suggest 
that we are seeing active star formation in galaxies from the upper end of the luminosity 
function. The apparent space density of Ha emitters in our survey is similar to that of 
the bright Lyman break objects identified at z > 3 and is within a factor of ~ 2 of that 
of present day L* galaxies. The grisms on NICMOS provide a first look at the enormous 
potential of low resolution near-IR spectroscopy from space. 
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TABLE 1 

The NICMOS Grism Survey Fields 
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J0923+8137 


32 


510 


1790 


J1326-4258 


19 


510 


1790 


J1958+4037 


05 


4080 


10740 


J0325-3618 


56 


1020 


3580 


J0923+8149 


32 


510 


1790 


J1326-4245 


19 


510 


1790 


J2000+2236 


03 


510 


1790 


J0349+1712 


28 


1530 


1790 


J0927-0502 


31 


1530 


5370 


J1336+5150 


63 


1530 


3580 


J2044-3124 


36 


2040 


5370 


J0354+0945 


32 


1020 


3580 


J0931-0450 


32 


510 


1790 


J1336-2959 


31 


1020 


3580 


J2044+6758 


15 


255 


895 


J0400-4433 


48 


1020 


3580 


J0931-0449 


32 


510 


1790 


J1337+7023 


46 


2040 


1790 


J2044+6758 


15 


1275 


895 


J0408-1207 


41 


510 


1790 


J0950+7222 


38 


510 


1790 


J1338+7009 


46 


510 


1790 


J2056-0443 


29 


1530 


1790 


J0408-1209 


41 


510 


1790 


J0954+6905 


40 


1318 


2556 


J1341+2824 


78 


4080 


14320 


J2106-2411 


39 


2040 


3580 


J0421+1926 


21 


1530 


1790 


J0955+6908 


40 


1510 


1917 


J1343+0218 


62 


510 


1790 


J2144+2813 


18 


510 


1790 


J0424+2707 


15 


510 


1790 


J0956+6937 


40 


1530 


3580 


J1343+0124 


61 


1020 


3580 


J2150+2845 


19 


1530 


5370 


J0425+1709 


21 


510 


1790 


J1020+3654 


56 


1020 


3580 


J1343+0225 


62 


510 


1790 


J2151+2854 


19 


1020 


3580 


J0430+2449 


15 


510 


1790 


J1039+4313 


58 


510 


1790 


J1345+0201 


61 


1020 


3580 


J2155+6334 


07 


1530 


1790 


J0431+1803 


20 


1530 


3580 


J1039+4145 


59 


6380 


3580 


J1346+0206 


61 


510 


1790 


J2155+6341 


07 


1530 


3580 


J0432+2424 


15 


1020 


1790 


J1044-5936 


00 


1530 


1790 


J1347+0219 


61 


1020 


3580 


J2155+6343 


07 


1530 


3580 


J0433+2526 


15 


510 


1790 


J1044-5948 


00 


3060 


7160 


J1402+5419 


59 


510 


1790 


J2156+6345 


07 


2040 


1790 


J0443+3612 


06 


1530 


3580 


J1053-0342 


48 


1530 


5370 


J1415+1910 


69 


2040 


5370 


J2157+6331 


06 


1318 


2556 


J0455+3029 


08 


510 


1790 


J1055-0418 


48 


1020 


1790 


J1432+0423 


56 


1530 


5370 


J2157+6338 


07 


7140 


5370 


J0458-0404 


26 


510 


1790 


J1056-0337 


48 


1020 


3580 


J1454+1218 


57 


510 


1790 


J2217-0336 


46 


510 


1790 


J0506+5250 


07 


2040 


3580 


J1057-0318 


48 


2040 


5370 


J1504+0110 


48 


510 


1790 


J2240+8031 


19 


510 


1790 


J0526-6602 


33 


2550 


3580 


J1120+2327 


69 


510 


1790 


J1507+0234 


49 


1530 


5370 


J2257+2050 
J2319+0809 


34 

48 


510 
1020 


1790 
3580 
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TABLE 2 
Emission- LiLNE Objects from the NICMOS Grism Survey 



Field 


Obj 




q(2000) 


^(2000) 


A 


^Ha 


flux" 


N„ 


W<,t, 


Wf;„ 


















^' 




erg/s/cm^ 




A 


A 


J0040+8505 


a 


00 


40 


43.93 


85 


05 


49.4 


1.727 


1.63 


1.5 


4.8 


250 


62 


J0055+8518 


a 


00 


55 


37.85 


85 


18 


17.3 


1.160 





76 


5.3 


5.2 


240 


25 


J0240-0140 


a 


02 


40 


10.85 


-01 


40 


53.0 


1.839 




80 


0.7 


12 


160 


49 


J0613+4752 


a 


06 


12 


59.93 


47 


52 


43.3 


1.341 




04 


1.3 


5.8 


60 


60 


J0622-0018 


a 


06 


22 


13.78 


-00 


18 


25.5 


1.402 




14 


2.9 


4.0 


403 


55 


J0627-6512 


a 


06 


27 


32.20 


-65 


12 


35.0 


1.742 




65 


1.8 


3.5 


176 


86 


J0738+0507 


a 


07 


38 


57.86 


05 


07 


35.2 


1.210 




84 


16 


10.9 


161 


32 




b 


07 


38 


58.37 


05 


07 


56.7 


1.77 




69 


0.9 


3.0 


180 


105 


J0741+6515 


a 


07 


41 


40.89 


65 


15 


23.8 


1.612 




45 


1.0 


7.4 


92 


72 




b 


07 


41 


45.37 


65 


15 


36.2 


1.879 




86 


1.8 


12.8 


165 


118 




c 


07 


41 


44.67 


65 


15 


06.9 


1.33: 




06 


1.8 


4.1 


879 


101 


J0917+8142 


a 


09 


17 


47.25 


81 


42 


21.4 


1.880 




86 


2.9 


9.4 


375 


131 




b 


09 


17 


44.96 


81 


42 


40.6 


1.883 




87 


0.9 


4.4 


224 


87 




c 


09 


17 


39.64 


81 


43 


06.5 


1.364 




07 


56 


115 


737 


32 


J0923+8149 


a 


09 


23 


34.32 


81 


50 


11.2 


1.572 




39 


4.2 


7.0 


547 


114 


J0931-0449 


a 


09 


31 


25.97 


-04 


49 


45.4 


1.299 




98 


24 


14.4 


550 


49 


J1039+4145 


a 


10 


39 


41.79 


41 


45 


04.4 


1.636 




49 


1.8 


10.7 


410 


118 


J1056-0337 


a 


10 


56 


56.41 


-03 


37 


03.9 


1.126 




72 


4.6 


3.1 


168 


36 


J1120+1300 


a 




20 


37.36 


13 


00 


07.1 


1.376 




10 


0.2 


6.8 


2809 


60 


J1120+2323 


a 




20 


24.29 


23 


22 


40.8 


1.557 




37 


1.9 


4.8 


660 


74 




b 




20 


22.36 


23 


23 


13.0 


1.397 




13 


1.4 


4.8 


456 


83 


J1134+0406 


a 




34 


04.71 


04 


06 


33.7 


1.260 




92 


5.4 


3.2 


1165 


86 


J1143-8036 


a 




43 


19.56 


-80 


36 


46.2 


1.538 




35 


1.2 


5.5 


310 


80 




b 




43 


19.73 


-80 


36 


46.9 


1.549 




36 


0.7 


5.4 


330 


53 




c 




43 


10.31 


-80 


37 


06.9 


1.668 




54 


1.2 


3.8 


353 


73 




d 




43 


18.46 


-80 


36 


30.2 


1.49 




27 


1.3 


4.1 


377 


298 


J1145+4733 


a 




45 


29.14 


47 


33 


35.2 


1.257: 




92 


1.1 


2.5 


250 


184 


J1218+4714 


a 




18 


36.42 


47 


14 


34.0 


1.56 




38 


0.7 


3.9 


108 


36 


J1237-I-6209 


a 




37 


49.04 


62 


09 


50.0 


1.802 




74 


0.9 


5.2 


172 


51 


J1237+6215 


a 




37 


38.85 


62 


15 


24.8 


1.778 




71 


0.4 


3.2 


108 


24 


J1237+6219 


a 




37 


06.57 


62 


19 


37.1 


1.557 




37 


1.8 


6.0 


607 


65 




b 




37 


05.28 


62 


19 


38.5 


1.348 




05 


1.3 


2.8 


546 


164 




c 


12 


37 


07.71 


62 


19 


24.8 


1.734 




64 


0.4 


4.2 


580 


303 



^The Ha line flux is in unites of 10~ ergs/s/cm . 

^Wiim is the 2(T observed equivalent width limit, which reflects the signal-to-noise ratio of the contin- 
uum near the wavelength of Ha line. 



-31- 



TABLE 3 
Photometry of Emission-Line Objects 



Field 


Obj 


Zfla 


H 


J-H 


logio(LHa) 


SFR 








mag 


mag 


erg s 1 


Mg/yr 


J0040+8505 


a 


1.63 


20.4 




42.40 


22 


J0055+8518 


a 


0.76 


19.7 




42.20 


14 


J0240-0140 


a 


1.80 


20.8 


1.1 


42.16 


13 


J0613+4752 


a 


1.04 


18.1 




41.90 


7 


J0622-0018 


a 


1.14 


19.1 




42.34 


19 


J0627-6512 


a 


1.65 


20.4 


0.8 


42.49 


27 


J0738+0507 


a 


0.84 


17.9 




42.78 


54 




b 


1.69 


22.0 




42.21 


14 


J0741+6515 


a 


1.45 


21.1 


0.8 


42.10 


11 




b 


1.86 


22.3 


0.7 


42.60 


36 




c 


1.06 


19.6 


1.0 


42.06 


10 


J0917+8142 


a 


1.86 


20.2 




42.81 


58 




b 


1.87 


20.2 




42.31 


18 




C 


1.07 






43.56 


324 


J0923+8149 


a 


1.39 


19.1 




42.69 


44 


J0931-0449 


a 


0.98 


19.0 


0.7 


43.11 


115 


J1039+4145 


a 


1.49 


20.2 


1.4 


42.39 


22 


J1056-0337 


a 


0.72 


18.8 


1.1 


42.10 


11 


J1120+1300 


a 


1.10 


20.8 


0.6 


41.14 


1 


J1120+2323 


a 


1.37 


20.5 




42.33 


19 




b 


1.13 


21.1 




42.01 


9 


J1134+0406 


a 


0.92 


19.9 


1.0 


42.40 


22 


J1143-8036 


a 


1.35 


21.4 




42.11 


12 




b 


1.36 


20.5 




41.89 


7 




c 


1.54 


20.8 




42.24 


16 




d 


1.27 


21.0 




42.09 


11 


J1145+4733 


a 


0.92 






41.71 


5 


J1218+4714 


a 


1.38 


20.3 


0.9 


41.90 


7 


J1237+6209 


a 


1.74 


23.3 


0.5 


42.24 


15 


J1237+6215 


a 


1.71 


20.8 


1.0 


41.87 


7 


J1237+6219 


a 


1.37 


20.5 




42.31 


18 




b 


1.05 


21.9 




41.91 


7 




c 


1.64 


22.3 




41.83 


6 



Note. — The Ha luminosities and star formation rates were 
calculated for Hq = 50km/s/Mpc and qo = 0.5. Here the 
listed star formation rates have not been corrected for the 
contamination from [Nn]6584,6548 line. As decribed in de- 
tail in the text, the correction to the star formation rates is 
29%. We only have J— H colors for some of the fields where 
the FllOW observations were available. "..." indicates no 
measurements. 
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TABLE 4 

Emission- Line Surveys in the Near-IR 



SITE 


Type 


^A 


Area 

n' 


Volume 
Mpc^ 


z range 


ergs/s/cm^ 


ergs/s 


W 


Ref 


Remarks 


Calar Alto # 1 


Imaging 


1% 


276 


30,000 


1.5- 2.8 


3.4 


3.8 - 19 


1 


1 


Targeted toward AGN 


Calar Alto # 2 


Imaging 


1% 


228 


20,600 


0.9- 2.4 


2.4 


1.0 - 9.5 


18 


2 


QSO metal line absorbers 


Keck/NIRC 


Imaging 


1% 


12 


815 


2.3- 2.5 


0.7 


2.5 -3.2 


13 


3 


QSO emission and absorbers 


ESO near-IR 


Imaging 


1% 


48 


10,000 


2.2- 2.3 


1.9 


6.2 - 6.8 


2 


4 


Radio Galaxies 


Lick/UCLA-camera 


Imaging 


1% 


180 


4255 


0.95-0.995 


3.9 


1.8-3.3 


2 


5 


R.Q. QSO emission z 


NICMOS Parallel 


GRISM 


45% 


58 


100,000 


0.8- 1.9 


0.4 


0.1- 0.9 


33 


6 


Random fields 



^All the flux limits are 3<t and in unites of 10""^^ ergs/s/cm^. 

'^The luminosity limits represent the minimum line luminosities which the survyes are sensitive to and are in unites of lO^^ergs/s. 

'^Number of emission line galaxies found in the survey. 

Note. — Co-moving volumes are calculated for Hq = 50 km/s/Mpc and qo = 0.5. 

References.— (1) Thompson et al. 1996; (2) Mannucci et al. 1998; (3) Teplitz et al. 1998; (4) van der Werf et al. 1997; (5) Teplitz et 
al. 1999; (6) This paper 



